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A GENERAL INTEGRAL FORM OF THE BOUNDARY-LAYER EQUATION FOR INCOMPRESSIBLE 
FLOW WITH AN APPLICATION TO THE CALCULATION OF THE SEPARATION 
POINT OF TURBULENT BOUNDARY LAYERS 1 

Bv Neal Tetebyin - and Chia Chiao Lin 


SUMMARY 

A general integral form of the houndary-layer equation is 
derived from the Prandtl partial-differential houndary-layer 
■equation. The general integral equation , valid for either 
laminar or turbulent incompressible boundary-layer flow, con- 
tains the Von Kdrman momentum equation, the kinetic-energy 
equation, and the Loifsianskii equation as special eases. 

In an attempt to obtain a practical method for the calculation 
of the development of the turbulent boundary layer, use is made 
of the experimental finding that all the velocity profiles of the 
turbulent boundary layer form essentially a single-parameter 
family. The general equation is thereby changed to a simpler 
one from which an equation for the space rate of change of the 
shape parameter of the turbulent boundary layer can be obtained. 

The resulting equation for the space rate of change of the 
■velocity-profile parameter is restricted by the assumption that 
the velocity profiles of the turbulent boundary layer can be 
approximated by power profiles. Two of the resulting equa- 
tions are used to calculate the distribution of the profile shape 
parameter over an airfoil for one experimentally determined 
pressure distribution. Although different assumptions were 
tried for the shearing stress across the boundary layer, the calcu- 
lated distribution of the profile shape parameter did not agree 
exactly with the experimental- distribution. 

An examination is made of the effect of using the experi- 
mentally determined single-parameter family of velocity pro-files 
■instead of the power profiles on certain functions that occur in 
the equation for the space rate of change of the velocity-profile 
parameter. One calculation of the distribution of the profile 
shape parameter over an airfoil is also made for the experi- 
mentally determined pressure distribution by using the single- 
parameter family of velocity profiles found from experiment. 
A comparison of the results with those of a calculation made 
with the same assumptions except for the use of power profiles 
shows some difference near the separation point. It is believed, 
however, that the apparent lack of reliability of the specific 
equations used to make the calculations is caused mainly by 
the lack of precise knowledge concerning the surface shear and 
the distribution of the shearing stress across the turbulent 
boundary layer. The present analysis emphasizes the need 
for information concerning the shearing stresses in turbulent 
boundary layers. 


INTRODUCTION 

An outstanding problem in aerodynamic theory is to cal- 
culate whether the flow will separate from the surface of a 
specific body and, if so, where the separation will occur. 
The concept of the boundary layer and the equations that 
describe the flow in it, introduced by Prandtl (reference 1) 
and first worked out in some detail by Blasius (reference 2), 
reduce the problem to solving the Prandtl boundary-layer 
equation when the flow is laminar. Because of the mathe- 
matical difficulty of solving the equation, approximate 
methods were developed for the calculation of the properties 
of the laminar boundary layer (reference 3). In some of 
these methods, for example, the Pohlhausen method (refer- 
ence 3) and the Wieghardt method (reference 4), a functional 
form is chosen for the velocity distribution through the 
boundary layer and is combined with either the Yon K&rm&n 
momentum equation alone (reference 5) or with both the 
Von K&rm&n momentum equation and the kinetic-energy 
equation (reference 4). The result is the replacement of the 
Prandtl partial-differential equation by one ordinary differ- 
ential equation in the Pohlhausen method and by two ordi- 
nary differential equations in the Wieghardt method. A 
solution of the ordinary differential equation or equations 
provides the boundary-layer velocity profiles along the body. 
These and other approximate methods that use. only the 
Yon K&rm&n momentum equation, or the momentum and 
kinetic-energy equations together, do not satisfy exactly the 
Prandtl boundary-layer equation. 

Because the flow in the boundary layer is more often 
turbulent than laminar in cases encountered in engineering, 
the problem of calculating the separation point is of even 
more importance for turbulent than for laminar boundary 
layers. In spite of the importance of the problem, however, 
less progress has been made in the development of methods 
for the calculation of the behavior of turbulent boundary 
layers than for laminar boundary layers. The lack of 
progress stems from the absence of an explicit independent 
equation for the shearing stress that is accurate enough to 
lead to a description of the flow when used with the Prandtl 
equation. 

The main attempts to obtain methods for the calculation 
of the behavior of the incompressible turbulent boundary 


1 Supersedes NAG A TN 2158, “A General Integral Form of the Boundary-Layer Equation for Incompressible Flow with an Application to the Calculation of the Separation Point of Tur- 
bulent Boundary Layers" by Neal Teterrin and Chia Chiao Lin, 1950. 
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layer in the presence of pressure gradients are those of 
references 6 to 12. The results of these attempts are unsatis- 
factoiy either because the assumptions upon which they rest 
are incorrect or because the equations used to make cal- 
culations were not derived from the boundary-layer equations. 

The analysis of reference 6 is based on the assumption that 
the velocity profile is a single-valued function of the ratio of 
the pressur’e gradient to the skin friction, an assumption 
shown to be incorrect by later investigators (for example, see 
reference 12). In the analyses of references 7, 10, 11, and 
12 the momentum equation is used, together with an auxil- 
iary equation, to calculate the distribution of velocity pro- 
files over a surface. In each of these four methods the 
auxiliary equation is not derived from the boundary-layer 
equations but is empirical. 

In reference 8, the equation that gives the variation of 
the mixing length across a pipe (reference 3) was used to 
calculate the_ velocity profiles. The fact that the mixing- 
length distribution across the boundary layer is not the 
same as across pipes is shown in references 13 to 15. 

Reference 9 does not provide a method for the calculation 
of the distribution of turbulent velocity profiles along a 
surface. It does, however, suggest that separation of the 
turbulent boundary layer always occurs when the numerical 
value of the nondimensional pressure gradient reaches an 
empirical constant. 

The purpose of the present investigation is to begin with 
the boundary-layer equation for incompressible flow and 
to proceed as closely to a method for the calculation of the 
behavior of the turbulent boundary layer as the present 
knowledge of the turbulent boundary layer permits. 

At first it might appear that the use of empirical auxiliary 
equations in methods for the calculation of the behavior of 
turbulent boundary layers can be avoided by developing 
a method similar to the Pohlhausen method which requires 
the solution of only the Von Kdrmdn momentum equation. 
For turbulent flow, however, in contrast with laminar flow, 
the conditions on the behavior of the velocity profile at the 
surface that can be obtained from the boundary-layer equa- 
tion seem to be of little or no use for the determination of the 
shape of the velocity profile across the boundary layer. 
Tins difference between laminar and turbulent flow makes 
inapplicable the Pohlhausen process in which a type of 
function is chosen to represent the velocity profiles, the 
function for the velocity profiles is combined with the Von 
Kdrm6n momentum equation, and the resulting ordinary 
differential equation for the space rate of change of the 
profile shape parameter is solved. 

An auxiliary equation for the calculation of the behavior 
of the turbulent boundary layer can, however, be obtained 
from the boundary-layer equation by making use of the 


experimentally verifiable fact (references 7, 10, 11, 14, and 

15) that all velocity profiles of the turbulent boundary 
layer form essentially a single-parameter family of curves. 
In the present analysis the Loitsianskii equation (reference 

16) is generalized by multiplying the Prandtl boundary- 
layer equation not only by an arbitrary power of the velocity 
in the boundary layer but also by an arbitrary power of the 
distance from the surface. The resulting equation is then 
integrated across the boundary layer and provides a general 
integral form of the boundary-layer equation, valid for 
either laminar or turbulent flow. This general integral 
form of the boundary-layer equation reduces to the Loitsi- 
anskii equation when the distance from the surface is raised 
to the zeroeth power, to the Von Kdrmfin momentum equa- 
tion when both the distance from the surface and the velocity 
are raised to the zeroeth power, and to the kinetic-energy 
equation when the distance from the surface is raised- to the 
zeroeth power and the velocity is raised to the first power. 

When use is made of the assumption of a single-parameter 
family of velocity profiles, the general integral form of the 
boundary-layer equation becomes a general equation for 
the rate of change along the surface of the velocity-profile 
shape parameter. This equation for the rate of change of 
the velocity-profile shape parameter is the desired auxiliary 
equation. 

The assumption of the single-parameter family of velocity 
profiles changes the problem from one of finding a solution 
of a partial-differential equation, the Prandtl boundary- 
layer equation, to one of finding a solution of two simul- 
taneous ordinarj^ differential equations, the equation for 
the rate of change of the shape parameter and the Von 
Kfirm&n momentum equation. The differential equation 
for the rate of change of the shape parameter, however, 
cannot result in a solution of the problem in the present 
analysis because a knowledge of the shearing stress is lacking. 
In the present analysis various assumptions are made for 
the distribution of shearing stress tlirough the boundary 
layer, and the distribution of the shape parameter over the 
surface of an airfoil is then calculated. Because of the arbi- 
trary assumptions for the shear distribution and the use of 
a flat-plate skin-friction formula, precise agreement between 
the calculated and experimentally obtained distributions 
of the shape parameter is not obtained. 

The problem of finding the shearing stress in the turbulent 
boundary layer remains. It is believed, however, that, if 
suitable approximations are found for the shear and surfaco 
friction, the equations presented herein should enable the 
development of the turbulent boundary laj’er to be calcu- 
lated with an accuracy sufficient for engineering purposes. 

The present work was begun while Dr. Lin was temporarily at 
the Langley Laboratory and was cont inued by correspondence. 



GENERAL INTEGRAL FOBS! OF THE BOUND ART-LATER EQUATION 


1069 


A 

a, b, c 

B 

v, 

H 

r 

g =T a 

gi 

gi 

=■4 

H 0 


SYMBOLS 

arbitrary positive integer in. shear polyn mial 

ri Qff 

coefficients in polynomial for I d£ 

exponent in expression for shear 
reference chord 
functional notation 


derivative of shear polynomial for \=0 
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exponent of y in derivation of general equation 

coefficient of u in equation for 8 

exponent in equation for power profiles f= f ” 
static pressure 
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radius of body of revolution 

coefficient of <f> in equation for 8 

velocity parallel to surface and at outer edge 
of boundary layer 

velocity parallel to surface and inside bound- 
ary layer, positive in direction of positive x 
velocity perpendicular to surface and inside 
boundary layer, positive in direction of 
positive y 
value of v at y—0 

coordinate parallel to surface, positive in 
direction from leading to trailing edge 
coordinate perpendicular to surface, positive 
outward from surface 

smallest value of y for which the difference 
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momentum thickness ^ /(I —j) dy ^ 


viscosity 



density 
shearing stress 
surface shearing stress 


ANALYSIS 


DERIVATION OF GENERA!. EQUATION 


The general equation is derived for the body oi revolution 
because the equation for two-dimensional flow can be 
obtained from this equation by letting the radius of a trans- 
verse section of the body of revolution become infinite. 
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The boundary-layer equation of motion for the body of 
revolution, also valid for two-dimensional flow, from refer- 
ence 3 is 


U 


bu . bu 1 dpi . 1 dr 

da: ' by p dx - ' p by 


(1) 


After multiplying through by u m , making use of the equation 
of continuity that is valid in the boundary layer of a body 
of revolution (reference 3) 

du . « dr a n 

bx^by r 0 dx~ 


and noting that 

n TT———~Sl 

p dx dx 

equation (1) becomes 

u m+1 / bu . by . u dr 0 \ . u da” 4-1 v du* 41 
m+1 V.da:"' by^~r 0 dx/^m + l bx "'m-fl by 



After equation (2) is written in a form in which each term 
vanishes at the outer edge of the boundary layer, each term 
of the equation is multiplied by y n and integrated from y— 0 
to y=8. The result (see appendix A for detailed develop- 
ment) is ; . _ 

L(m + 1 )]+£ § [N-n(J-M)]-nQ $ 

b— 

— (3) 


Equation (3) is the general integral form of the boundary- 
layer equation. 

The Von K4rmdn momentum equation is obtained from 
equation (3) by letting m=0 .and ft=0, the equation for 
kinetic energy is obtained by letting m=l and n=0, and 
the equation for moment of momentum is obtained by 
letting m=0 and n=l. 

In the case 77t=n=0, 


Nd=J‘(l-f)fdy=e 


or 

A r =l 

Also, 

r s 


Le = (J-l)dy 


Jo 

or 

1 

II 

‘ S,|® 
1 

II 


n approaches 0; thus 

nQ8 n =nJ y n ~ 1 (l—f Bt+1 )dy 

=[V(1 ~/ m+1 >]*- jV ^ (1 —f m+1 )dy 

The first term drops out if n^O. Then, by taking the limit 
n— >0, the result is 

lim 71(20*= 1 


n-)0 


Hence, when m=n=0, equation (3) becomes 

dd . 8 dU i 0 drb ro 

dx^U dx ^ ±Jl) *r 0 dx U~jLP 


(4> 


Equation (4) is the momentum equation for flow over a 
body of revolution with flow through the surface. For two- 
dimensional flow, equation (4) becomes 


d 0 , 0 dXJr-rr I ftp r q 


(5> 


when the value of ^ from equation (4) is substituted into 
equation (3), the result is 

L _ NH }-ljg n [(J-AI)+W-$[nQ-N(n+l)] 

== “^[ (m+1) r^ m,7,, |f rf,7+Cri+1)iNr ] (6) 


where 


r 



The assumptions contained in equation (6) are the usual 
boundary-layer assumptions. Equation (6) is valid for 
both laminar and turbulent flow. 


FORM OF EQUATION (6) FOR SINGLE-PARAMETER FAMILY OF 
VELOCITY PROFILES 

Equation (6) is now to be placed in a form valid when the 
velocity profiles form a single-parameter family of curves 
(/=/(» !,£?)). For this purpose the term h of equation (0) 
is modified in the following manner: 

By definition, 

Because? depends only on y and 11, 

d(l -f) b(l-f) by b(l-f)dH 
bx by bx~*~ bH dx 


It can be easily verified that all the integrals, except Q, 
involved in equation (3) have finite integrands as n ap- 
proaches 0. The limit nQ, however, approaches unity as 


From the definition of y 


by _ y dd 
bx~ 8 dx 
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then 


Jo &* V eixj o’ b, » + ia:Jo dff V 

or, after an integration l>y parts of the first term on the right- 
hand side, the result is 

dH HbJ , 
dxjo m dy 

The expression for then becomes 


But 


and 


^r»- a ~ fmn) ^’M iy ) iy 

fV(l —f m+1 )dy=M8* +1 

fV (1 —f tn+i )fdy=N6 n+l 


Then with 


the present no restrictive assumptions regarding the shear 
are made. The form of the kinetic-energy equation for a 
- single-parameter family of velocity profiles is obtained from 
equation (8) hy placing m=l and n=0 and dividing by 
N'(H)=K'(B). Thus, 


dH K(H- 1) 6 dU | 

<K+2 \ 


1 T 0 

dx" K! U dx 

\ K' ) 

pC/ 2 


'K—l\v 0 
. K' )U 


(9) 


The symbol K represents the ratio of the kinetic-energy 

thickness (1 — j f)f iy to the momentum thickness. Note 

that in the derivation of equation (9) from equation (8), the 
assumption of a single-parameter family of curves is not 

restricted to the case-^=0. 


RESTRICTION OF GENERAL EQUATION TO POWER PROFILES 

The data in figure 1 show that the power profiles defined 
by /=£* are a good approximation to the “standard” profiles 
derived by fairing experimental data (reference 10). Equa- 
tion (8) can be further developed by using the assumption 
that /=f p . After some fairly lengthy calculations (see 
appendix B), equation (8) becomes 

B — 4j?(p + l)(2j>+l)[p(m+2)+n,+ lI 8_ dU 

- dx pm+n- f-1 U dx~^ 


the expression for I x can be written as 

< 7 > 

When the expression for I x from equation (7) is substituted 
into equation (6) and equation (4) for ^ is used, the follow- 
ing equation is obtained: 

s fr(S?+" 7 Mf i [- n (H+ i )y-ifl+K«+ i )+ 

N(E—m)\ M)-N- 

(m+l)J"V”^§J^]+ 

^[n(J-M)-N+nQ] (8) 


where has been used. Equation (8) for 6 

is applicable both to two-dimensional flow and to flow over 
a body of revolution. 

In equation (8) all the integrals, except the one involving 
the shear ratio g, are functions of H, m, and n only. For 


[p(m+2)+n+l]£ J-IJ-+ 

2p(p+l)Lp(TO+2)+n-H] d 0 

p(m+l)+n 77 ^ ' 

As a first approximation the assumption has been made that 
j—l v even when 

The occurrence of the arbitrary positive integers m and n 
in equations (8) and (10) requires an explanation. In order 
to determine why m and n appear, equation (8) is written in 
a different form. By making use of the definitions for A 7 , 1, 
J, M, L, and Q and integrating by parts where necessary in 
order to eliminate terms that contain ij* -1 , the result is 


g+»/=-(m+i)j"W- ( 4 |-| i. 


-n(H+ 1 ) ( J- M) + L)m -f- 1 )■ + N(H—m) 

=(m+l) [(ET+1) §£ -Jjdv+f-ljdv 


JT-tm+l) J’Tf’V ^ in 
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n(J-M)-A^+nQ=(m + l)J o S/ \ tt / m ^ (l- J # 7<*u) dr, 


Equation (8) then becomes _ _ 


dx 




( 11 ) 


By iiaing the assumption of a single-parameter family of curves directly in the. partial-differential equation (1), the 
following ordinary differential equation is obtained: 


( 12 ) 


-°-l 
n dH P U* [ 

@fj 

O’+S 0+uf 



,®o 

'U 

H -> s=- 
zO I/O 

1 1 

[V^-O] 

6 dx ~ 






The concept of a single-parameter family of velocity pro- 
files is consistent with equation (1) and with particular 
functions for tq/pU 2 , g, and / when the right-hand side of 
equation ( 12 ) is independent of 17. When the right-hand side 
of equation (12) is independent of r,, the right-hand side of 
equation (11) is independent of m and n. Equations ( 11 ) 
and (12) are then identical. 

To obtain an equation fox 9 that does not contain either 

m or n or both, the functions r 0 IpU 2 , g, and / must there- 
fore be such that the right-hand side of equation (12) is 

independent of r , ; the solution of the equation for 6 then 

provides a solution of equation (1). Note that the problem 
is to find a solution not of equation (1) alone but of equation 
(1) and the independent relation for the shearing stress in 
turbulent boundary layers; this relation is at present un- 
known. 

The nature of the approximation made in the^gresent 
analysis, in order to obtain a specific equation for may 

dH 

be clarified by noting that a specific equation for is ob- 
tained from equation (8) by choosing the functions T 0 /pZP, 
g, and / and substituting an arbitrary positive integer for m 
and an arbitrary positive integer for n. The calculated 
distribution of Iff over a body for arbitrarily chosen functions 
for TolpU 3 , g, and /is then consistent with the momentum 

equation and one of the integral equations for For 

example, if m = l and n=0, both the. momentum and the 
kinetic-energy equations are satisfied but no other ones. 
If m=0 and n= 1, only the momentum and the moment of 
momentum equations are satisfied. In the present analysis 
only the momentum, the kinetic-energy, and the moment of 
momentum equations — equations which have familiar phys- 
ical meaning — are used. 

As noted previously, equation (11) is independent of m 
and n if the functions rofpU 2 , g, and/ are such that the right- 
hand side of equation (12) is independent of 17. In this case 
a solution of equation (1) results and the func tions.ro/pC/ 3 , 
g, and / and the calculated distribution of H satisfy every 
particular equation obtainable from equation (11), (10), or 
(8) by assigning positive integers to m and n. 

Note that m and n cannot both be made zero in equation 
dN 

(8) because 0 for m=n=0. .. If m and n are both 



n 

(a) JT- 1.310, 1.518, and 1.710. 

(b) H=1S80, 2.057, and 2.277. 

(e) H-2.455 and 2.701. 

Fioube 1 .— Comparison ol eiperimental and power velocity proflloa. 
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zero, equation (8) becomes 0 =e= 0. It is also noted that 
equations (8) and (10) are valid both for flow over a body of 
revolution and for two-dimensional flow. 

For m=l and n— 0, equation (10) leads to the equation 
for kinetic energy 

»ir- 

(H+ S -S=l£r | «) $,+ 

(H+l) (3S - — 1) c n 
4 U 

where the relation for power profiles 2p+l=H has been 
introduced. This form of the energy equation can also be 
obtained from equation (9) by noting that, from the defini- 
tion of K and the equation for power profiles, 

K 2(2g+l.) 4 S 

3p + l 3S-1 

A comparison of the values of K obtained from this formula 
and obtained from the standard profiles is given in figure 2. 

The equation of moment of momentum for power profiles 
is obtained from equation (10) by letting m — 0 and n= 1; 
it is 

n dH g(ff+l)(g*-l) 0 dl! , 

° dx 2 U dx~' 

(fl ra -l)[l?+(H+l)J o I r|f:dr]^+(S s -l)^ (14) 

In this equation the term involving the shear distribution 
may be rewritten as follows: 

It then involves the mean shear inside the boundary layer. 

ATTEMPTS TO DEHIVE A RELATION GOVERNING THE CHANGE OF THE 
FORM PARAMETER 

In most of the recent analyses of the development of a 
turbulent boundary layer, an empirical relation governing 
the change of the form parameter H is usually introduced. 



It is clear that equation (10) automatically furnishes such 
relations if the shear distribution is known. In this section, 
three attempts are described to establish such a relation, 
These attempts are based on the following simple assump- 
tions for the shear distribution: 

(a) The shear distribution depends only on — which 

is equal to the Pohlhausen parameter multiplied by a factor 
(reference 3). 

(b) The shear is constant across the boundary layer. 

(c) The shear distribution depends only on the form 
parameter of the velocity distribution. 

The first two assumptions are used either with the energy 
equation in the forms given by equations (9) and (13) or 
with equation (14) for the moment of momentum. The 
last assumption is used with equations (13) and (14) jointly. 

(a) Shear distribution depending only on the Pohlhausen 
parameter. — The first assumption follows the original idea 
of the method of Yon Kfirm&n and Pohlhausen in using 
polynomial approximations together with the boundary 
conditions obtained by successive differentiation of the 
equations of motion (reference 3). Fediaevsky (reference 8) 
appears to have been the first to introduce it into the inves- 
tigation of turbulent boundary layers. When the shear- 
stress through the turbulent boundary layer is assumed to 

be a polynomial of fifth degree inf =y satisfying the follow- 
ing boundary conditions: 
at y~0 


at y=S 


the following expression is obtained: 

^=d-r) 3 [i+(3+x)r+3(2+x)r 2 i (is) 

The shear distribution g is a function of f and X, where 

^ 5 dpi 6 dU 5 pU 2 

— t 0 dx ‘ U dx ~d r 0 


dr dpi 

T ~ r ° dy dx 

O 

II 

/O I/O 

r=0 0 

by 

*1=0 

by* 


The particular boundary conditions at y = 0 restrict this 
development to the case-^=0. 

From the shear distribution (equation (15)) the calculation 
may be made of the coefficients P and S. The attempt to 
calculate P and S' by using the standard profiles together 
with equations (9) and (15) was, however, unsuccessful for 
two reasons. First, the ratio 5/0, which must be known, 
could not be accurately determined from the standard pro- 
files. Second, for reasonable values of 5/0, the calculated 
values of P were positive for values of H for which P should 
be negative. 

The calculation of the part of P independent of the shear 
profile was then made both for the standard profiles and the 
power profiles by making use of the kinetic-energy equations 
(equations (9) and (13), respectively); the comparison is 
shown in figure 3. The closeness of the results suggests 
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that it is permissible to use power profiles as an approximation 
for calculating P and S. From equations (13) and (15), 

griwggwB+w ] < 16 > 

and 

S— (8H— i)[g— ( a4 ( ^ 7 7g +9) ] dri 


The functions P and S, given by equations (16) and (17), 
respectively, are shown in figure 4. 

The fact that the equation 

•§-r*+s* 


where P and S are obtained from equations (16) and (17), 
respectively, does not predict the behavior of the turbulent 
boundary layer is shown as follows: _Xet w=0; then, for 

H greater than approximately 1.5) should be negative. 

Because S from equation (17) is positive for £T> 1.2, it 

follows that is positive. This conclusion is incorrect; 
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Fioori 4.— Values of P from equation (15) and values of S from equation (17). 


therefore, the function for S (equation (17)) is inconsistent 
with the known behavior of turbulent boundary layers. 

To show that the function for P (equation (16)) is in- 
consistent with the known behavior of turbulent boundary 


layers, let Ha* 1.4. By making positive and largo, 
becomes positive and large because P given by equation 


(16) is positive. 


dU 


For positive values of -g— 


however, it is 


known that should be negative. The function for P 

(equation (16)) is therefore inconsistent with the known 
behavior of turbulent boundary layers. 

In order to determine whether functions for P and S that 
do not result in obviously incorrect conclusions can bo 
obtained bj' making the shear polynomial satisfy a greater 
number of boundary conditions at the outer edge of the 
boundary laj’er, the shear polynomial is generalized by 
writing 


g=( i - ir [i ■ + Ai+ x f (i - sr a + At) 

(18) 


The boundary conditions at the surface that are satisfied by 


equation (18) are 



. 5=1 

or 

r=r 0 

II 

or 

dr_dpi 
dy~ dx 

^=0 

dt 2 

or 

*I»0 


At y=d, the conditions that are satisfied are 


5=0 


l g 

df A ~ 1_ 


0 


In order to evaluate the integral f f* ~ df in equation (13) 
r Jo Of 

*g. 

or 


. 

the term is written as 


where 


Of 


=5i + \?s 


5i=d -f/ [a+a(a + Dir ] -a (i-rr-> [i + An r s ] 

and, for A ^ 1, 


5 S =(l-f) 4 (l + 2Af)-fA(l + x4f)(l-fy 1 - 1 
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By using the expression for g u the equation obtained for /S' is 

S=2$p+iy- 


(2p+l 

Alp r i 

1 A ■ 

A(A+ 1)(2>+1) 1 ) 

Up+i 

(A+p)(A+p-l)(A+p-2) . . . (p+l)b 

1 A+l+iP 

2(A+ 2+p) ( A+ 1 +p)J) 
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(19) 


By using the expression for g 2 the equationfor P is found to be 

P=-2(2p+l)(3p+l) s - • • (P+1) [a+T +p + (A+2+p){A+ 1 +p)] j (20) 


To avoid positive values for S obtained from equation (19) 
for H<3, it is found that A must be 1 in the expression for 
gi. It is also found that, to avoid positive values for P in 
equation (20) for £T> 0, A must be “ in the expression for 
gt. The values for S and P then become 


O (3ff-l)(H- l)(H-3) 
H+5 

P=-H(H-1)(3H-1) 


( 21 ) 

( 22 ) 


The expression for P (equation (22)) is the same as the 
coefficient of jj in equation (13) ; letting A— > makes the 

coefficient of X in equation (18) become zero. The shear 
profile then contributes nothing to the coefficient of jj 
in equation (13). 

Equations (21) and (22) for S and P, respectively, were 
tested by making a computation of H and 0 for the pressure 
distribution given in table I of reference 10. The computa- 

tion began at -=0.075 with the vaLues given in table I of 
c 

reference 10. The equations used are 

e g=-g(g- i)(3H- i) w + (3 ^~ 1) -^ ^ )(g — } * 


H+5 


and 


dd 


£=-(H+ 2)c+* 


The equation for is the Von KArmdn momentum equation. 

The equation for <j> was obtained from reference 17 and is 

0.006535 


4>~- 


F+ 8 


The calculated distribution of H along x was far from the 
experimental curve. 

In an attempt to reduce the sensitivity of the equation for 

6 to the shear distribution, the moment of momentum 
ax 

equation (equation (14)), in which the shear appears in the 
coefficient of <f> only as a mean value, is used. "When the 
generalized expression (equation (18)) is used for the shear 
distribution g, the result obtained is 


'f=[- 


H(H— 1 )(H+ 1) 2 , 3H(H+ l) s 1 , 

2 ' t ~(A+2)(A+3)J c,!i ' 

CEP— i) * 


(23) 


where \{f = 0 as required by equation (18). To keep the 
coefficient of w negative for all positive values of H, A must 
equal =° in the coefficient of «. The shear distribution is 
then independent of the pressure gradient. To make the 
coefficient of <f> negative for values of H near 3, A must have 
the smallest value that it can take; therefore, let A=1 in 
the coefficient of <f>. Equation (23) then becomes 

n dH — H(EP— 1)(H+1) ■ (H 2 — 1)(H— 3) f , 

9 1Z = 2 4 * (24) 

A calculation for the example in table I of reference 10 with 
equation (24) resulted in a computed curve for H that was 
far from the experimental curve. 

(b) Assumption of constant shear across the boundary 
layer. — All the computations of H have led to values of H 
much larger than the experimental values. Therefore, in 
order to reduce the calculated values of H it is necessary to 
increase S. In order to increase S, the assumption of con- 
stant shear across the boundary layer is made. For constant 
shear it can be shown that 


/: 




by letting g=(l — T) s and taking the limit of the integral 
as B~* 0. Equation (14), after the assumption of constant 
shear is' introduced, becomes 

u -(pp- l)*+CH*-l)* 

In order to make ^jj?=0 at H=1.286 for the coef- 

ficient of <t> was arbitrarily changed to H 2 — 1.286 2 . The 
equation then becomes 


m = -H{HAr\)(IP 1 ) M _( g _ 1-65 3 fr+(H»-l)^ (25) 
ax & 

This equation was used for the computation of H with 
^=0, and the results’ for the example given in table I of 
reference 10 are shown in figure 5. 

The assumption of constant shear across the boundary 
layer was also combined with the kinetic-energy equation. 
When the power profiles and the assumption of constant 
shear are used in equation (13), the kinetic-energy equation 
becomes 


9^=-H(H-l)(3H-l)c 
(H+1)(3H— 1)_, 


(H— 1)(3H— 1) 
2 


£+ 
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Figure 6 .— Comparison of eiperimental values of H and values calculated by equation ( 2 S). 
f. JiBfaH) . s~m- U658; »-0. 


The function — 2?(2?— 1)(327— 1) is shown in figure 3 and 
the function — — — — — >in figure 6. When the stand- 
ard profiles are substituted for the power profiles and the 
assumption of constant shear is made, the kinetic-energy 
equation (equation (9)) becomes '7' 

n dH K(H- 1) K ~ 2 / K - 1 , 


dx 


K' 


K! 


<t>- 


K' 


J6/JT 1 ") 

where the function -g, — - is shown in figure 3 and the 

2 .... .. ... . 

function — - is shown in figure 6. The results of these 


calculations of 27 (with ^=0) are shown in figure 7. In this 
case, the use of power profiles makes the result somewhat 
different from that obtained by using the standard profiles. 

(c) Determination of S by the simultaneous use of the 
energy and moment of momentum equations. — It seems 
obvious that, if equations (13) and, (14) were exact, the 
coefficients of «, <j>, and ^ in equation (13) would be equal to 
the coefficients of as, <f>, and \(/ in equation (14). The ratio of 
the coefficients of « is 
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F;gub& 6.— Values of S from kinetic-energy equation with the assumption of constant shear. 

For power profiles, for eiperimental profiles, =77-. 

2 it 



Figure 7 —Comparison of experimental values of II and values of II calculated by kinetic 
energy equation with constant shear across boundary layer. 

The ratio of the coefficients of $ is 


27(27+ l)(H*-l) 

2 ' ‘ (27+1) 2 

-H(H — 1 ) (327- 1 )‘ “ 2(323 — 1 ) 

(H+l ) 2 

The curve of 2 ( 3 %— 1) S* ven “■ % ure 8 and is seen to be 
dose to unity. 


4(iT 2 -l) 4(2?- 1) 

(H+l)(8flr-l)“ 327-1 

The curve of is also given in figure 8. The values 

are far from unity for small values of 2? but become equal 
to unity for 2?= 3. 
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Figure 8. — Ratios of pressure and suction, (actors In energy and moment of momentum 

equations. 

Equating the coefficients of <j> results in 

<B*-1) [b+(B+ 1)JV if] 

or 

(ff-xKH+i^V | | it 

=B ! (3 — B) (26) 

Now let 

J\^dt=a+bH+cH* 

and 

When the shear distribution is assumed to depend only on 
B, the integrals in equation (26) are functions of B alone. 
Because equation (26) is then an identity, the coefficients 
of the various powers of B can be equated to zero. The 
resulting equations are: 


For B> 


for H 1 


for IP 



— a— b + 3j— 1=0 


a— b— c— |i+3g— ^-=3 


a+b— c — g 2+33=— 1 


forB 4 


6+c-|3=0 


for IP 


The results obtained, are: 


c=0 


a= 


6 = 


28 

128 

180 


128 


c=0 


56 


2 — 


T28 

32 


Therefore, 


and 


X 


o or 


128 

40 

128 

7 — 4517 


32 


X' 


„ 2 ?>g 


dt- 


7+4B+5B 2 


is 

Equation (14), the moment of momentum equation, becomes 


6 


a > — 


dH B(B*-1)(B+1) 
dx 2 

(B— 1)(3B— 1)(7+22B+15B 2 ) 
32 


*+(ff-l)* (27) 


and equation (13), the energy equation, becomes 

dg=-ff(H-l)(3S-l) a - 

(B-1)(3B-1)(7+22B+15B*) , , (B+ 1)(3B-1) , 
_____ *+ * 

(28) 

The variation of 27 with a: for the initial values and the pres- 
sure distribution given in table I of reference 10 was computed 
by using a modified form of equation (28). In order that 

g=0 at values of B in agreement with experiment when 

&>=0, the coefficient of ^ in equation (28) was replaced by 


(B-B 0 ) (3B-B„) (7+22B+1 5 B 2 ) 
32 


where 


Bo=B 0 (i?«) 

The variation of B 0 with if# was calculated from the equation 
log i0 B 0 = 0.5990— 0.1980 log 10 i? s -0.01S9 (logic if#) 2 
which was derived to represent a faired curve through the 


for IP 
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experimental data (see fig. 9); the data were obtained from 
reference 13 and from British results that are no b generally 
available. The result of a computation of H for \p=0 and 
with equation (28) modified as follows 

0^=-H(ff-l)(3H-l>- 

(H-H <i )(2H-H 0 )a+22E+15H i ) , , (H+l)(3ff-l) , 

32 ~<f>+ 4 

(29) 

is given in figure 10. 

Assumptions (b) and (c) lead to somewhat better results 
than assumption (a) although they ar.e still not as satisfac- 
tory as those obtained from the purely empirical relations 
introduced in references 10 and 12. ...It is clear that this 
difference is caused partly by the inaccuracy of the simple 
assumptions about the shear distribution and can be im- 
proved by using better descriptions. However, in view of 
the limited present knowledge of the shear distribution, it 
does not seem worth while to make more complicated 
assumptions. 



Figure 9.— Variation of H% with R#. 


It. may be noted that the final equations obtained for the 
change of the form parameter by the three assumptions are 
all of the form 

e^-=P{H)*+S{H)4> 

where m=-f? w- and d >=— ^ This form is used in refcr- 
U ax pU * 

ence 12, but a different form is used in reference 10. 

INVESTIGATION OF ENEBGY EQUATION 

Since none 'of the three assumptions for the shear dislribu- 

dH 

tion results in a dependable equation for an investigation 
is made to determine whether a result common to the three 
assumptions — namely that the coefficient of.-^jjin the equa- 
tions for fr- is a function of II alone — is very far from truo 

by using experimental data and the kinetic-energy equation 
without any assumption for the shear. 

If no assumptions other than the boundary-layer assump- 
tions are made and if in equation (6) n = 0 and m= 1, the 
result is 


(30) 

If the assumption of a single-parameter family of curves is 
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made (J=J^,B)), then K=K(B), and equation (30) becomes | or 
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or, for ^=0, 


K+2 


C U, M_ 

Jo 


in 


K' 


-a 7~y/ K + f- 


K' 


K' 


e~=<s>k(jS){Z+&) 


(31) 


where 


nm- 


-(H-l)g 


K' 




If the assumption is made that g=g(y,3), then /=/(g r 3) 
and J f dg is a function of 3 only. Therefore, £q=£o(22). 


Equation (31) then becomes 

e^=<t>k(H)[it+um\ 

In order to obtain an estimate of the quantity 1 


(32) 


2 r i 


under the assumption that f=f(g, 3), reference 12 is used. 
Equation (7) of reference 12 may be written as 


JH 


8-^=4>ki(H) [£— 2.065(22'— 1.4)] 


(33) 


where 


k^3)=d 


— .SCff-l.* 


Note that Gamer’s equation (equation (33)) has the form 
the kinetic-energy equation takes when the assumptions 
that and that g=g{ v, 22) are used in the kinetic- 

energy equation. The kinetic-energy equation (equation 
(31)) can also be placed in the form of equation (32) when 

the more general assumption that F Q (ij,H)+F 1 (ri,H ) 

is made for the shear distribution. For the purpose of 

2 f* 1 

obtaining an estimate of the value of 1 tJ. / dg, the 

quantity £-4-£o(22) in equation (32) is assumed to be identical 
with the quantity £ — 2.065(22— 1.4) in equation (33). 
Then 

£o(H)= -2.065(22-1.4) 
and for 27=1.5, for example, 

£ 0 (22) = —0.2065 

therefore, 

o ( 1 - iJT / *) _ “ 0 " 20 ' 5 


1—IJV dg= — 0.1032 

Therefore 1— -j^J / dg is the difference between two quanti- 
ties, each of which is much larger than their difference. It 
follows that, in order to determine £<>(22) for values of 3 not 
close to separation with any accuracy, } and g must be 
known with relatively good accuracy. 

It may be noted that the moment of momentum equation 
is also sensitive to g. This sensitivity can be seen by writing 
the coefficient of <f> in equation (14) as 

(EP-1) [(ff+1) (l- JW)-l] 

When it is noted that the integral f g d$ is of the order of 

unity and that 3 lies between 1.2 and 2.6, the sensitivity of 
the coefficient of <j> to g becomes clear. 

In an attempt to determine whether £ 0 is determined 
mainly by 3, all the data that were used in reference- 10 were 
used to compute £ 0 by making use of equation (32) in the 
form 

•m 

The surface-friction coefficient <j> was calculated by the 
formula (from reference 17) 


4>= 


0.006535 


2fu I/B 


and k (22) was calculated by the expression obtained from the 
moment of momentum equation 


*(fl)= 


ffCEP-DCff+l) 

2 


The values of £ 0 plotted against 3 are given in figure 11. 
The effort to determine whether £ 0 is a function mainly of 3 
is inconclusive. At least part of the scatter occurs because 

(p? and ^r~ were obtained from curves faired through ex- 
ax ax 

perimental points. In addition, the calculation of £ 0 requires 

e d3 
8 dx 

the subtraction of £ from aEL operation which further 

decreases the accuracy of the calculated values of £>. 

DISCUSSION 

Although equation (6) is valid whenever the boundary- 
layer assumptions are valid, the equations for that 
result after additional assumptions are made do not lead to 
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Figure 11,— Values of (t plotted against H. 

good agreement with experiment. The first of the additional 
assumptions made is that all velocity profiles of the turbulent 
boundary layer belong to a single-parameter family of curves. 
The experimental data of references 7, 10, 11, 14, and 15 
substantiate this assumption. 

The second assumption is that the single-parameter family 
of curves can be approximated by power profiles. The data 
in figure 1, in which velocity profiles are compared, and also 
the data in figures 2 and 3, in which K and P are compared, 
show this assumption to be good, at least for i?<1.8. 

From the data in figures 1 to 3, it is inferred that power 
profiles can be substituted for the standard velocity profiles 
without greatly affecting the calculated distribution of H 
against x for H<C1.8. To. test' this inference, the kinetic- 
energy equation was used with the assumption of constant 
shear across the boundary layer; the result is shown in 
figure 7. As expected from the data of figures 1 to 3, the 
effect of the substitution of power profiles for the standard 
profiles is noticeable only for £T>1.8. It thus appears that 



the inaccuracy of the equations for 


that were 


caused mainly by the surface-friction law that was used 
and by the assumed shear distributions rather than by the use 
of the power profiles. 


The data of references 12 and 15 show skin frictions that 
increase, strongly in the region upstream of the separation 
point before dropping to zero at the separation point. On 
the other hand, the skin-friction data presented in reference 
14 indicate that the skin friction falls monotonically to zero 
as the separation point is reached. In the present analysis 
a skin-friction law obtained from experiments on flat plates 
is used. It is therefore probable that part of tho inaccuracy 
in the equations used to calculate H is caused by the use of a 
relation for the skin friction that does no t give correct values 
when there are pressure gradients along the surface. 

The assumptions for the shear distribution that were made 

to obtain a specific equation for were 


(a) The shear distribution depends only on the ratio of the 

, 5 dpi & ca 

pressure gradient to the skin friction — or— — 


(b) The shear is constant across the boundary layer 

(c) The shear distribution depends only on the form param- 
eter of the velocity distribution 

Because none of these simple assumptions is derived from a 
knowledge of the details of the turbulent flow, it is not likely 
that any of them are valid. When it is recalled that the 
coefficient of £ in both the kinetic-energy and the moment of 
momentum equations is sensitive to the shear distribution, 


it is not surprising that a reliable equation for 


Ur IX J 

was not 


found. 


In order to obtain a reliable equation for 


from equation 


(8) it thus seems necessary to calculate the surface shear and 
the shear distribution across the boundary layer more accu- 
rately than in the present analysis. Efforts should therefore 
be made to understand the mechanics of turbulent shear 
flow sufficiently well to provide an independent relation for 
the shearing stress that will predict the behavior of turbulent 
boundary layers when used with the Prandtl boundary -layer 
equation (equation (1)). 


CONCLUDING REMARKS 


A general integral form of the boundary-layer equation is 
derived from the Prandtl partial-differential boundary-layer 
equation. The general integral equation, valid for either 
laminar or turbulent incompressible boundary-layer flow, 
contains the Yon Ktirmdn momentum equation, the kinetic- 
energy equation, and the Loitsianskii equation as special 
cases. 

In anattempt to obtain a practical method for the calcula- 
tion of the development of the turbulent boundary layer, use 
is made of the experimental finding that all tlio velocity 
profiles of the turbulent boundary layer form essentially a 
single-parameter family. The general equation is thereby 
changed to a simpler one from which an equation for the space 
rate of change of the shape parameter of the turbulent 
boundary layer can be obtained. 




GENERAL INTEGRAL FORM OF THE BOUNDARY-LAYER EQUATION 


1081 


The resulting equation for the space rate of change of the 
. velocity-profile parameter is restricted by the assumption 
that the velocity profiles of the turbulent boundary layer can 
be approximated by power profiles. Two of the resulting 
equations are used to calculate the distribution of the 
profile shape parameter over an airfoil for one experimentally 
determined pressure distribution. Although ' different as- 
sumptions were tried for the shearing stress across the bound- 
ary layer, the calculated distribution of the profile shape 
parameter did not agree exactly with the experimental 
distribution. 

An examination is made of the effect of using the experi- 
mentally determined single-parameter family of velocity 
profiles instead of the power profiles on certain functions 
that Occur in the equation for the space rate of change of the 
velocity-profile parameter. One calculation of the distribu- 
tion of the profile shape parameter over an airfoil is also made 


for the experimentally determined pressure distribution by 
using the single-parameter family of velocity profiles found 
from experiment. A comparison of the results with those of 
a calculation made with the same assumptions except for the 
use of power profiles shows some difference near the separa- 
tion point. It is believed, however, that the apparent lack 
of reliability of the specific equations used to make the calcu- 
lations is caused mainly by the lack of precise knowledge 
concerning the surface shear and the 'distribution of the 
shearing stress across the turbulent boundary layer. The 
present analysis emphasizes the need for information concern- 
ing the shearing stresses in turbulent boundary layers. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Ya., May 22, 1950. 


APPENDIX A 

DETAILED DEVELOPMENT OF EQUATION (3) 
Equation (2) can be written so that terms of the form 


v, mJrl — U m+l appear explicitly; therefore, .each term will 
vanish at the outer edge ' of the boundary layer. The 
resulting equatjon is 


1 buU” +l , 1 bvU m+l , 


+ 


m + l dx m- f-1 dy 

1 (u m+1 —U m+1 )u dr Q , 1 uU m ~ l dr ft 


m + l 


•J-o 


dx m+l r 0 dx 


mT r dl . U m dr 

<=u n b ~~j — 1 5 — 

dx p oy 


(Al) 


or, after simplification, 


Tl (Z7“ +1 — u* +l ) 0 ]— 


m+l 

1 ({7 m+I — u m+l )u dr 0 


m + l 


r 0 


dx 


•Cu Si U—uU m ) ^-=— — 
{u 1 dx p by 


(A2) 


Equation (A2) is now multiplied through by y a and inte- 
grated with respect to'y from y= 0 to y=S. The resulting 
equation is 


fa* ^ (U^-u^u dy- 
TlJo ¥ * 5y ' tt m+1 )p dy— 

< j^j o y*(u m U—uU m ) u m y* ^ dy 


m 


m 


or, after simplification and substitution of the formula for 
the differentiation of a definite integral 

J* y n ^ (U m+1 — u x+! )udy=-j^j y n (JJ m+1 —u m+T )udy 

the following equation results: 

Tl Jo ^ T/ (U n+1 —u m+l ) v dy— 

r*- 

[(#)"+] y" ‘tvjfyr % d « 

(A3) 

By integration by parts, 

J* y n ~ (U m+1 — u m+r )v dy=— nj (U a+1 — u m+l )vy n ~ 1 dy 

and equation (A3) becomes 

- 5 TFT E p " + ’JV[ 1 -(r)‘ +, ]r‘ i >'+ 

JL-j\u^i- U «+r) vy *-i dy- 


m 

m 
U m+ 1 


1 U m+1 dr 0 n 


M-csna*- 


m + l r 0 dx, 

v ’*'^§l\& m ^ y ' d y=)$yy'% iy (A4) 


218837— S3 69 
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The velocity v can be eliminated from the term 

f (Cl”* 1 — u m+1 )vy n ~ 1 dy 
Jo 

by the following development: 

The velocity v may be written as 

v= !’U dy+ ''‘ 

or, by use of the equation of continuity, 
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Now let . _ 


• b(U-u) J dU , 1 dr a C B , , U dr 0 
Si-^-Vis+jr.TSi, <P-u>iy--^v+„ 

r > \ r jv~u) ,dU /■, u\ , il 7. 

“J. tF a* l “3x t 1 ~v)\ dv ~ y 

U dr 0 rv/ u\, U dr, ~ 
r o da; Jo \ II) y r, dx 

^fe°Jo ( 1 ~v) dy + Vo 

or, after terms are collected and / is substituted for u/U, 
the result is 

The term (U m+1 —u m+1 ) vy*~ x dy can now be written as 
j* (U m+1 —u m+t )vy n ~ 1 dy— J (f7" +l — 'ii m+1 )y* -1 

! u jo dy+ (j^ + r 0 ii)\_ fo {1 -a d v~y ]+ v » j d v 


or 


r 

Jo 


(U m+1 — u n+l )vy*- 1 dy 


- c ,-«/;r-a-y-«)[/;^V 5 ^]rfy+(f+ 


773?) ^ + '/>-‘H-?’ + ')[/ t '(i-3)3 s -!/] J,+ 

o,U m+1 £ y n ~ l ( 1 —j m+x )dy 


j>*-a -/-■)[/; 

JV-*(i -/”+’) [J o ' (1 -J)dy\ dy=J6" +1 

r y* (1 -y* +1 ) dy=Md n+1 
Jo 


and 


£ y n -\l - f n+1 ) dy= Q8 n 
The term J (U n+1 —u m+1 )vy*~ 1 dy now becomes 
J (27 s 1+1 — u m+1 )vy’ l ~ l dy=U m+2 I 1 d n+1 -{- 


Now let 


and 


£y n {l-f M+x )fdy=N8^ x , 
£ US*) y*dy=Ld*+ x 


Equation (A4) can then be written as 

-5VFT lx I 

n / dU , U dr. 


771 + 1 


(w+Ei) C7 " + ‘ y - M#,+,+ s4-i^ £7 - +, « s '- 

~* + * j^NB n+1 + 27”* 1 Ld*+ X =— dy 

m + 1 j,. dx dx p J 0 J a by y 

(AS) 

c? 

After ^ 27* +2 iV0' ,+1 is expanded and terms are collected, 
equation (A5) becomes 

( „ +1)Jv g +e (f -»/,)+ . 

^ g [N(tti +2) -n (J —M)—L(m+ 1)] + * 




kli 1 w-»w-ioi-ne^-(»+o^J^>v- 5 a<i, 

where 


= y. 


7? = 



APPENDIX B 

SIMPLIFICATION OF TERMS IN EQUATION (8) FOR POWER PROFILES 


CALCULATION OF ^+nl 
act 


■The definition of N8 n+l is 

iW+i- r (1 dy 

s-r" [° 


therefore, 

and 


or 

The definition of l8 nJrl is 


therefore, 


Then 


Id* +l =£(l—f n+l ) ^ dy) y n ~ 1 dy 


n£\ 1 -f.^(J’IL d7 ,y- iJ , (BD 

By integration by parts, 

• n j 0 ^ d n)v n ~ 1 d7i^ 

- X‘" [ C1 M- <m + X” • <**> 

When equation (B2) is substituted into equation (Bl), the 
following equation is obtained: 

Use is now made of the power-law assumption 




[hen 


a/ i-2 P 2 , r log r 

. dH~2(p+l)(2p + l) i ^ 2 


and 


d v J 0 bH av J bH 


l-2g» 


M 


2(p+l) 2 (2p + l) 


ir ^Mogj; ^ -R 

r + 2L p+1 (p+l)‘J 


.[2(2?+ 


—2tf 


2(p+l)(2p+l) 


f s, 4 


r 2j log t 


dZV. 


] 


After a lengthy manipulation, is found to be 

diV . r m+1 (p+l)' t (2p+l)"(pm-l-7i) 

2 p* +1 [ 2 (m+ 2 )+n+l\ l 

where use has been made of the following equation: 

S (g + l)(2p + I) 
d p 

CALCULATION OF _V 

The definition of Nd a+1 is 


Nd a + l =£ Cl -/ m+I )/r dy 


When /= £ ” is used, the equation for N is 

p(m+l) 


W4 V+1 


or 


d) [p(m+2)+?i,+ll(p+7i.+l) 

• (m-H)(p + l)- i - 1 (2p+l) n+1 

2?“[p(m-f-2)+7i,-t- 1] Cp+»-t 1) 

CALCULATION OF J 

From the definition of J0 n+1 , 

Je s+1 = J (1 — /® +l ) (1 — f)dyj y n ~ 1 dy 

When/= the equation for J is 

'-GTr,*-*™* 1 ( r_ l?i) t "' dt 

or, after a lengthy manipulation, 

j. fr+l)*(2?+l)**‘ ( » +I) 

P* A 

f (p+l+7i)[pm+n+2(p4-l)]+g(m-H)+l+7t ) 

l(7i+ 1)I?+ l+»)[ 2 Km + I)+n+l][p(TO+l)+p+l+»]) 
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CALCULATION OP SI 

From the definition of A/0’ 1+1 

A/0" +1 = [ (i — f m+v )y n dy 
Jo 

When/= X 9 , the equation for M is 

M=(m+1) 0>+» w (3g+l>* +' 

CALCULATION OF L 

From the definition of Z0 8+1 

Zfl ,,+1 = f (1 —f*- l )fy*dy 
Jo 

When /= X* , the equation for L is 
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By substitution and simplification 


L=(m— 1) 


Q>+l)‘ +1 (2g+l) ,,+1 ‘ 

p n (p+n+l)(pm+7i+l) 


CALCULATION OF V(g — > 

dh +nI 

From the expressions for J and M, the expression for 
J— Mis 

j-jtf- SElim+ir* 1 (w+1) 

tp + l+m)[pm+n+2(p+l)] + [p(m+l)+l+n] 
(n+i)(p+l+n)[pCm+l)+«4-l][p(m+l)+p+l+n] 

| n (P+l) a+1 (2p + l) a+1 ■ 

' p" (n+1) [p(m + l)+7i+l] 

After a lengthy simplification, the result is 

t ~(m + l)(p+I)»(2p+l)^ 1 . 

p a (p+l+7i)[p(m-f'i)+p+l+n] 
or 

r_ Tfwg | 2n(m + l) (p+ l) n+1 X2p+ l) ,i+1 
' ^ p* (p+ 1 +n)\p(m-\- l)+p + 1 +n] 

where N=2p-\-l was used. The expression obtained for 
N(H— m) + (m+l)Z is 

N{H— m)+(m+l)Z 

(m + l)(p + l) lt+1 (2p+l) l,+1 2p(pm+m+7i) 
p s (p + n+l) [p (m+2)+n+l](pm+»+l) 

and the expression obtained for 

N(H-m)-n(J-M) (H+ 1) +lm+l)L 
is 

N(H—m) —n(J— A/) (H+ 1) + (m + 1 )L 

_ 2(m + l)(g+l)’ t + 1 (2p + l) n+1 [ft(pm+n+l)+p(pm+CT+7i)] 
P"(p+ l+n)[p(m + 2)+n+l](pm+7i+l) 


N(H — w) — n(«7 — A/)(Z/+ l}+(m+l)Z 

dN, : 

w +nI 

_ -4p(p + l)(2p + l)j>(m+2)+n+lJ 
pm + n + 1 

It can also be shown that 

— A/+n(J— A/)— (m + 1) f ri n f n ^~dy 
Jo oy 


dN, T 
dH +nI 


2 P 


pm+n 


[p(m + 2) + n + l]|2p+l + tp(m+2)+n+l]^ 


f. if ) 


EVALUATION OF -A’+n U-.UJ+nQ 

From tlie results for N and J— M 

y | w. . -im + l)(p+l)»(2p+l)»^ . 

A+n(J A/1- p » [pim + 2)+n+l] 

For (?) 'the development is 

Q8«= Va-r^dy 

Jo 

and with.f=X v > the following expression is obtained for np^Q: 
m + I (p+l)'*(2p+ll y 


Q- 


n p”- l [p(m + l)+n] 


Then, by substitution and simplification, for n?^0, 

p"Lp(m + 2)+n+l][p(m + l)+n] 

(B») 

If use is made of the previously derived result that nQ= l 
for n— 0, the following equation is obtained for n= 0: 


-N+nV-U)+«Q- f £$l\ 


(B4) 


If n is placed equal to zero in equation (BR), equation tB4) 
results; therefore, equation (B3) is valid for n— 0 as well as 

715^0. 

Then, for all values of n, 

— N+n(J— Ai)+ng _ 2(p+l)p[p(m+2)+n+l] 
dN, T p(m+l)+n 

dH +nl 


GENERAL INTEGRAL FORM OF THTD BOUNDARY-LATER EQUATION 
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